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We study the low-temperature properties of linear Josephson-j unction arrays capacitively coupled 
to a proximate two-dimensional diffusive metal. Using bosonization techniques, we derive an effective 
model for the array and obtain its critical properties and phases at T = using a renormalization 
group analysis and a variational approach. While static screening effects given by the presence 
of the metal can be absorbed in a renormalization of the parameters of the array, backscattering 
originated in the dynamically screened Coulomb interaction produces a non-trivial stabilization 
of the insulating groundstate and can drive a superconductor-insulator transition. We study the 
consequences for the transport properties in the low-temperature regime. In particular, we calculate 
the resisitivity as a function of the temperature and the parameters of the array, and obtain clear 
signatures of a superconductor-insulator transition that could be observed in experiments. 
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I. INTRODUCTION 

Low-dimensional superconductors are systems display- 
ing a surprisingly complex and rich physics, allowing the 
study of paradigmatical phenomena in condensed mat- 
ter physics, like quantum phase transitions and quantum 
critical behavior, electronic localization. Coulomb block- 
ade, etc.^^ In particular, an intriguing superconductor- 
insulator phase transition (SIT) was observed experi- 
mentally in superconducting film^^M^, wire^^'f^, and in 
ultrasmall-capacitance Josephson junction arrays (JJAs) 
in two^ and one dimensions'^ giving rise to an in- 
tense theoretical activityP^^ In this transition, as one of 
the parameters is varied (e.g., the normal-state resistance 
of the film, the thickness of the wire, the Josephson cou- 
pling Ej in the array, etc.) the groundstate of the system 
changes from superconducting to insulator. 

In one-dimensional (ID) superconductors, one partic- 
ular kind of excitation, the so-called quantum-phase slip 
(QPS) processes, have been recently the focus of an in- 
tense research.^ The interest is based both on the puta- 
tive role of QPS in the SIT in ID'^, as well as for their 
potential uses in novel qubit architectures'^, a fact that 
has stimulated recent interesting experimental research 
in IDJJAs ;^^'^'' A phase-slip is a discrete process occur- 
ing in a ID superconductor, in which the amplitude of 
the order parameter vanishes temporarily at a particu- 
lar point, allowing the phase of the order parameter to 
change abruptly in units of 27r. In particular, a QPS is a 
phase-slip excitation originated in macroscopic quantum 
tunneling of the phase of the order parameter. 

On the other hand, since the seminal works by Caldeira 
and Leggett"^, it has been known that dissipation in 
macroscopic quantum systems plays a central role. For 
instance, in a two-dimensional JJA capacitively coupled 
to a proximate two-dimensional electron gas (2DEG), 
Rimberg et al. observed a tunable- SIT upon variation 
of the backgate voltage Vg applied to the 2DEG.'^ In 
that work, it was shown that Vg has the effect of tuning 



the sheet resistance in the 2DEG through the mod- 
ulation of its electronic density, a fact that in turn mod- 
ifies the electromagnetic environment of the JJA. It was 
argued later by Wagenblast et al?^ that due to the in- 
complete screening of the Coulomb interaction provided 
by the 2DEG, the 2D plasma mode in the array is over- 
damped and the charging energy in the junction Ec is 
renormalized to higher values, producing a SIT when- 
ever the ratio Ej/Ec ^ 1, with Ec the renormalized 
Eq' While this scenario is reasonable in a 2D geom- 
etry, in a IDJJA capacitively coupled to a 2DEG, the 
screening provided by the metal is typically very efficient, 
and a significant damping of the ID propagating plasma 
mode is not expected. This leads to the naive conclu- 
sion that in the 1DJJA/2DEG geometry, a dissipation- 
driven SIT should not occur. However, in a closely re- 
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Figure 1: Scheme of the system under study. The IDJJA is 
capacitively coupled to the metallic film, which controls the 
electromagnetic environment. A gate voltage allows to modify 
the sheet-resistance of the film, modifying the dissipation 
in the IDJJA. 

lated Luttinger-liquid system placed in proximity to a 
metallic plane, a dissipation-driven quantum phase tran- 
sition was predicted to occur. This transition is driven 
by backscattering events originated in the Luttinger liq- 
uid under the effect of the dynamically screened Coulomb 
interaction. It is therefore interesting to study to what 
extent the same dissipative processes will affect the dy- 
namics of QPS in ID superconductors in proximity of a 
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diffusive metallic plane. Indeed, the question of dissipa- 
tion in ID superconductors is an active area of research, 
and theoretical predictions point towards the important 
role of intrinsic and/or extrinsic dissipatio n me chanisms 
in determining their T = phase diagramP^I^ 

In this article we explore the possibility of a 
dissipation-driven SIT in a IDJJA capacitively coupled 
to a proximate diffusive 2DEG. We concentrate in par- 
ticular on the low-temperature phase diagram and on 
the transport properties of the array. Using a bosoniza- 
tion approach, we derive the dissipative effective action 
from a microscopic Hamiltonian, and we elucidate the 
role of dissipation in the SIT. One important conclusion 
in our work is that for weak dissipation, the transition oc- 
curs always between a superconducting and an insulating 
phase, in contrast to oth er works predicting quadrupolar 
and normal phases .1^^^ We believe this is a consequence 
of a different kind of dissipation in the model. We find 
that, except in the experimentally challenging situation 
in which the Cooper-pair density in the superconducting 
islands is not commensurate, the SIT is al ways of the 
Berezinskii-Kosterlitz-Thouless (BKT) type,^^^^ and is 
originated in the unbinding of QPS/anti-QPS pairs. Dis- 
sipation stabilizes the insulating groundstate through the 
introduction of friction in the dynamics of the ID super- 
fiuid density, a fact that could be observed experimen- 
tally in the dc-resistivity of the IDJJA. Specifically, we 
predict a resistivity of the form g (T) AiT^^ + 
in the superconducting phase, and q{T) ~ rfFe^l^ j 
in the insulating phase, with A the insulating gap and r] 
the dissipation parameter. 

The paper is organized as follows. In Section |TI| we de- 
rive the effective model for a IDJJA coupled to a 2DEG, 
in Section |III| we derive the T = phase-diagram as a 



function of the parameters of the model. Section IV 



devoted to the study of the experimental consequences 
of our results, and finally in Section |V| we present a sum- 
mary and our conclusions. 



II. MODEL 

We start the analysis by considering an ideally isolated 
J J A, with length L ^ 00. To simplify the analysis, we 
neglect in the following the fermionic degrees of freedom 
forming the Cooper-pairs at a microscopic level. This 
"boson-only" approximation is belived to describe cor- 
rectly the critical properties of a JJA at temperatures 
T <C Tc, with Tc the superconducting critical temper- 
ature in the bulk of the superconducting island.^^ The 
usual description of the isolated, infinite IDJJA is given 
in terms of the quantum phase model^ 
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order parameter Oi at every site i in the array. These 
variables obey the usual phase-number commutation re- 
lation in the BCS groundstate, i.e., [Oi^nj] = iSij.^'^ The 
first term in Eq. ([T]) represents the charging energy, with 
Vij the unscreened Coulomb interaction [cf. Eq. (|3|] 
between the excess charges at sites i and j, and n cor- 
responds to an average charge imposed, e.g., by external 
gate voltages. The second term is the Josephson energy 
contribution, parametrized by Ej. In the following we 
use the convention h = ks = l- 

The critical properties of model Eq. ([T]) are more con- 
veniently studied using a field-theoretical approach, valid 
for fluctuations of wavelengths much larger that the lat- 
tice parameter of the array We therefore introduce 
the coarse-grained superfluid density Sp{x)^ defined as 
Sp{xi) = {rii — n) /a, and we expand the Josephson term 
as Ej cos {Oi — Oj) :^ Ejo^ [VO{xi)]'^. At low tempera- 
tures, the continuum limit of Hamiltonian Eq. ([T]) reads 



^JJA ~ 2 J V [x — x' ^ 0) 8p (x') 

+ ]^Eja f dx {VO{x)f . 



(2) 



Here the ID superfiuid interacts via the hare Coulomb 
potential, which we define for convenience as 
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(3) 



(1) 

The dynamical variables of this model are the number 
of Cooper pairs rii and the phase of the superconducting 



where r = |r| and z are, respectively, the distance in the 
x?/-plane and along the azimuthal direction between two 
point-charges (cf. Fig. [T]). Here the lattice parameter a 
acts as the short-distance regularization of the interac- 
tion and er is the permitivity of the insulating medium 
surrounding the islands. Note that in Eq. ([2| we do not 
assume an a priori short-ranged, screened interaction as 
is usually done when dealing with JJAs.^^ This will re- 
sult as a natural consequence of the interaction with the 
2DEG (see below). One problem of this field-theoretical 
approach is that Mott-instabilities (crucial when the su- 
perfiuid density is commensurate with the lattice) are 
lost in Eq. ([2| after taking the continuum limit. One 
way to cure this problem is to introduce a phenomeno- 
logical term Hi = — J dx Vi (x) p (x), where the effective 
superfiuid density p {x) [cf. Eq. ([lO|] couples to the 
phenomenological potential Vi (x), having the same peri- 
odicity of the lattice. 

The electrons in the 2DEG are described by the Hamil- 
tonian 

(7 ^ 

^ ^ + i d^vd^v' Sp2B (r) V (r - r^ 0) Sp^B (r ) , (4) 

where the fermionic field-operator r]l = r]l (r) creates 
an electron in the 2DEG with spin projection a at spa- 
tial position r = (x,?/), and Fimp = "^mp (r) repre- 
sents a weak static impurity potential which provides 
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a finite resitivity and dissipation in the metal. In 
terms of 77J (r) , r]a- (r), tiie density-operator p2D (r) in tiie 
2DEG writes p2D (r) = Y.a (r) Va (r), and Sp2B (r) = 
p2B (r) — Po,2D, witii po,2D tiie average density in tiie 
metal. 

Finally, the interaction between the IDJJA and the 
2DEG placed at a distance d (cf. Fig. [T]) is described by 
the Hamiltonian 



^int = J dP^dx' Sp {x^) V {x^ - r, d) Sp2B (r) . 



(5) 



Our goal in this Section is to derive an effective model 
for the IDJJA capacitively coupled to the 2DEG. To that 
end we introduce the partition function of the system^^ 



where S is the Euclidean action of the problem 



S = 5'jJA + ^2D + ^int, 



(6) 



where 



n(3 n n(3 

SjjA = dr dx idrO {x, r) p (x, r) + / dr Hjja (t) , 
Jo J Jo 

S2D = J dr J d'^rf] (r, r) {dr - P2d) V (r, r) + H2D (r) 



Sint = dr Hint (r) . 



Here /i2D = k^/2m — eVg is the effective chemical po- 
tential in the metal, with = Ikpl the Fermi wavevec- 
tor, and Vg the gate voltage applied to the 2DEG, which 
allows to change the value of po,2D, and therefore, the 
sheet-resistance Ru- 

The first step in the derivation of an effective model 
for the array is to integrate out the fermionic degrees 
of freedom f] (r, r) , r] (r, r) in the 2DEG. Assuming that 
the term Si^t can be treated perturbatively (we check the 
consistency of this assumption later), the integration of 
the fermionic degrees of freedom in the metal yields 



eff :^ 5'jjA — ^ \ drdr' J dxdx' Sp (x, r) 

X Vscr {x -X ,T - r ) bp {x ,t') . (7) 



We do not provide the details of this derivation here, and 
we refer the interested reader to Refs. 1251 and [26l In Eq. 
^ we have introduced the ID effective screening poten- 
tial 'User [x — x' — t'\ which encodes all the screening 
effects provided by the 2DEG. This quantity writes more 
conveniently in Fourier representation aP^ 
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W2D (k,0) X0,2D (k,Wm)' 
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where uo^a = 2iTm/P are the bosonic Matsubara 
frequencies,^^ and k = {k^k±) is the wavevector in 
2D, where we have made explicit the component kj_ 
in the 2DEG, perpendicular to the IDJJA. The quan- 
tity V2D (k, d) = (27reVer) exp (- |k| VdPT~a^) / |k| is 
the 2D Fourier transform of the Coulomb potential Eq. 
([3|. We assume that the length of the array is L < ^loc, 
with ^loc the Anderson localization length in the 2DEG, 
a condition well fulfilled in practice. In that case, the 
density-density response function in the 2DEG, aver- 
aged over disorder configurations, writes Xo,2D (k, a;„ 



J), where D and ^f2D respec- 



tively, the diffusion constant and the density of states (at 
the Fermi energy) per spin projection.'^ 

We can now define the total effective retarded interac- 
tion 



^eff {k, OJm) = ^ID {k, 0) - VscT {k, UJm) ■ 



(9) 



where vm {k^O) = 2e'^Ko {\k\ a) /cr is the Fourier trans- 
form of Eq. ([3| in ID, and Kq {x) is the zeroth- 
order modified Bessel function.'^ Physically, the effec- 
tive potential Vef^ (/c, Um) describes the interaction among 
charges in the array, both via the instantaneous interac- 
tion viB (^,0) arising from the direct intrawire Coulomb 
interaction, as well as indirectly via the coupling to the 
diffusive modes in the 2DEG, which corresponds to the 
'retarded interaction v^cr {k^Um) Eq. (|8|. 

We now introduce a more convenient representation 
of the superfluid density in the IDJJA. To motivate our 
approach, we first note that in the absence of Josephson 
coupling [i.e., Ej = in Eq. ([T])], the Cooper-pair oc- 
cupation number is a good quantum number in each 
island, fixed by n via the application of an external gate- 
voltage. Increasing Ej will evidently introduce fluctua- 
tions in rii due to the transfer of Cooper-pairs between 
neighboring islands, and is no longer a good quan- 
tum number. However, we expect that in the experi- 
mentally interesting regime Ej/Eq ~ 1, where is the 
characteristic charging energy in the island, fluctuations 
Arii = rii — n will be of order An^ ±1, and that all 
other charging states such that |Ani| 1 will be en- 
ergetically forbiden. We therefore truncate those states 
from our description and focus on charge-fluctuations of 
Aui = ±1. In terms of a continuous field ^(x), which 
is slowly varying on the scale of a, the superfluid density 
in this effective model can be more conveniently written 
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(10) 



where the parameter po is defined as po = 1 /a in the com- 
mensurate case. Note that po is an effective parameter of 
our model, and cannot be interpreted as the total phys- 
ical density, in contrast to truly ID systems. Only the 
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fluctuations 5p{x) = p{x) — po have a physical meaning 
in our model. 

In order to obey the phase-number commutation rela- 
tions in the BCS-groundstate,^^ note that the field <p {x) 
must verify the new commutation relation 



[6> (x) , V0 (xO] =i7rS{x-x'). 



(11) 



The contribution in squared brackets in Eq. (10) de- 



scribe long- wavelength density fiuctuations around the 
average value po, while in the last term, each contribution 
describes low-energy density fiuctuations of momentum 
k ~ 2ppo^ where p is an integer. When replaced into the 
effective action Eq. ^ we obtain the following effective 
model 



Sef^ [(/>] = ^0 [^] + ^1 



(12) 



where the contribution 5*0 corresponds to a Luttinger liq- 
uid modeP^ 



^0 [^] = 
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resulting from the slow fiuctuations of the density 
Sp{x) ~ — (x) /tt and from the hydrodynamic (i.e., 
{k^ujm} 0) sector of '^eff (^, ^m)- u and K are respec- 
tively the velocity of the ID plasmon and the interaction 
Luttinger parameter 



K 



IT 



u = a^/EJEo, 



(14) 
(15) 



where Eq = e^/2Co is the charging energy with respect 
to the ground, with Co = e^a / Ahi {2d / a) the effective 
ground capacitance of the Josephson junction. In our 
treatment, due to the screening provided by the 2DEG, 
the static effective potential v^f^ (/c, 0) is effectively short- 
ranged for distances x ^ d^ and therefore the Luttinger 
parameter is a constant. ^^In terms of the capaci- 
tance matrix Cij of the IDJJA, this amounts to neglect- 
ing the interjunction capacitance C, since this contribu- 
tion, although relevant for density fiuctuations of mo- 
mentum k ^ a~^, drops off in the long- wavelength sec- 
tor /c ^ (i.e. the interaction is screened in a length 
Lgcr ~ ayC/Co).^In the present context, the Luttinger 
parameter K physically represents the competition be- 
tween coherence and charging effects in the array [cf. Eq. 
(p!4|)]. Therefore, a large parameter K favors supercon- 
ducting correlations, while a small value of K tends to 
destroy superconductivity due to strong charging effects.^ 
The dissipative parameter r] is defined as 



c Ru 
erSTT Rq '' 



(16) 



where R^ is the sheet-resistance of the 2D film and c is a 
numerical constant of order c ~ ^(1)- Eq- (13) with 



a non- vanishing 77 describes a ID plasmon- mode with 
a finite lifetime F ~ \k\ /vjW^ Physically, a broadening 
of the ID plasma mode occurs due to coupling to the 
diffusive modes in the 2DEG. The term ~ ri\uOm \ \k\ in 



Eq. (13) is the result of combining the leading contribu- 
tion in powers of \ujm \ / Dk^Y \ k\ (with /ctf the Thomas- 
Fermi momentum in the 2DEG) in the expansion of the 
retarded potential '^eff (^, ^m), and the long-wavelength 
fiuctuations of the density ~ (V^)^, which contributes 
a term ~ k'^ |0(A:,cj^)|^ in Fourier representation. Note 
that since the scaling dimension of the term ~ \k\ \uOm\ 
is 2, the critical properties of the system are not modi- 
fied. Moreover, for a metallic plane [cf. Ref. [10] with 
Ru ~ O.li^Q, T] 10~^ <C 1 and it can be effectively 
ignored, allowing us to write 
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\^{k,uj^)\\ (17) 



This assumption greatly simplifies the analysis, since the 
action Sq recovers Lorentz-invariance in space-time. 

The next term Si in Eq. (12) originates in the phe- 
nomenological potential Vi (x), which has the same peri- 
odicity of the array. Therefore, it can be decomposed in 
Fourier components as Vi {x) = Vn cos (Qnx) , with 
Q = 27r/a. In general, all terms other than p = n = 
in Eq. (10) are rapidly oscillating and vanish under the 
integral sign. However, if Qn = 27rppo^ or equivalently 
poa = n/p (i.e., the average density of bosons is commen- 
surate with the lattice), then the term ~ J dx Vi (x) p {x) 



yields a term e 



-i{Qn—27Tppo)x 



1 which is not oscillating 



and, in addition to the term p = n = 0, we have the 
additional term 

A 



^1 M 



aro 



/ 



dxdr cos {2(j){x^r)) ^ (18) 



where we have only kept the most important commensu- 
rability {p = 1), and where we have defined the dimen- 
sionless parameter A 



A = ViTo. 



(19) 



and the short-time cutoff tq = a/u. The term Vq can 
be reabsorbed in a redefinition of the chemical potential 
of the externally imposed charge n, so in the following 
we will not consider it. Physically, the dimensionless pa- 
rameter A is related to the QPS rate in the Josephson 
junction by Fqps = A/tq. Estimated experimental val- 
ues for Fqps are in the order of ~ 1 GHz^^, which yields 
A 0.06. 

The final term in Eq. (12) comes from the dissipa- 
tive part of Vscr {k^uom)' Due to the strongly oscillating 
factors ~ g-*27rppo^E Eq ^q)^ results in the local 



dissipative term 

S2[(t)] = -^ j dxdrdr' 



El cos 2p [(j) (x, r) 
p 

p>0 ^ 
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(20) 
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This contribution is consistent with that of Ref. ^7^, ob- 
tained in the context of Luttinger liquids capacitively 
coupled to diffusive metals. In spite of the small mag- 
nitude of 7^, we will show that this contribution has im- 
portant consequences for the critical properties of the 
IDJJA, in contrast to the term proportional to rj in Eq. 
(pl. 



In the following we study the critical properties and 
phases of the model obtained in Eq. (12). 
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III. PHASE DIAGRAM 

A. Weak-coupling renormalization group analysis 

We first focus on the phases of the IDJJA at T = 0. 
To that end, we perform a weak-coupling renormaliza- 
tion group (RG) analysis of the model Eq. (12), as- 
suming that Si and ^2 in Eqs. (18) and (20) respec- 



tively are weak perturbations to the Luttinger liquid So 
in Eq. (17). Since the action is Lorentz-invariant 



in space and imaginary time, we adopt an RG proce- 
dure that rescales homogenously space and time. As 
usual, we assumme that the original theory is defined 
up to a certain momentum cutoff A (/) = Age"^ (with 
Ao ~ and we study how the action 5*0 is renormal- 

ized upon integration of high-energy modes in a window 
between A (/) /s < | q| < A(/), with 5 = e"^ , 
have employed the compact notation q = {/c, 
X = {x, ur}. 

We obtain the perturbative RG-flow equations of the 
model by performing a one-loop correction in S2 and a 
two-loop correction in ^i, and requiring that the term 
5*0 is invariant upon scaling. We obtain the RG-flow 
equations 



where we 
-^} and 



dK{l) 

dl 
du{l) 

dl 
dX{l) 

dl 
dT]{l) 

dl 



-2Trri (l) - {2TTf K (l) X"^ (l) C K^l) , (21) 

-2Trri{l)u{l)K{l), (22) 

[2-K{l)]X{l), (23) 

[1-2K{1)]^{1), (24) 



where the numerical constant C is of order unity. 

Note that both and 5*2 tend to destroy supercon- 
ducting correlations in the Luttinger liquid phase, a fact 
that is reflected in Eq. ( [21] ) where the Luttinger parame- 
ter K (l) is renormalized to smaller values, meaning that 
charging effects are enhanced. This can be interpreted 
as an effective increase of the charging energy Eq in Eq. 
(14). In addition, since 5*2 is the only term that breaks 



Figure 2: Schematic phase diagram of the IDJJA in the K — rf 
plan e, obtained from the integration of the RG-flow Eqs. ( 21 )- 
(24), with the initial parameter Ao = A (/ = 0) = 0.01. An 
increase of Ru in the 2DEG, and consequently, of the dissipa- 
tive parameter 77, can induce a SIT. Note that, in absence of 
dissipation, the critical value Kc — Kc (Ao) — 2.1 is slightly 
shifted with respect to the value Kc (Ao ^ 0) ^ 2. 



When K (I) < 2, the perturbative parameter A (/) 
flows to strong-coupling [cf. Eq. (23)], and the per- 
turbative RG procedure is no longer valid. In the limit 
77 ^ we recover the usual Mott-transition of the BKT- 
type described by the sine-Gordon model, and below the 
critica l value K c = 2, the IDJJA is in the insulating 
phase. Using Eq. (14), this means that in ab- 

!|2l 



the Lorentz invariance of the theory, note that the plas- 
mon velocity u{l) is proportional only to r]{l)^ and is 
independent of A (/). 



sence of dissipation, the SIT occurs for Ej/Eq = (2/7r) 
Note that our situation corresponds strictly to the case 
when the superfluid density in the IDJJA is commen- 
surate to the lattice, and is in clear distinction to the 
non-commensurate situation (i.e., A = 0), where dissi- 
pation (i.e., the term S2) becomes relevant for K (1) < 
1/2, inducing a different kind of non-superconducting 
groundstate.^^ 

In the present case, the scaling dimension of the dissi- 
pative parameter rj (1) is always smaller than that of A (/), 
which means that for K (l) 2, is a stronger perturba- 
tion as compared to 5*2. Therefore, one would expect the 
nature of the non-superconducting groundstate to be de- 
termined essentially by ^i. However, based on this fact, 
one could naively conclude that the term S2 is unimpor- 
tant near the SIT, a conclusion we prove incorrect. In 
fact, a more detailed analysis reveals the importance of 
the term 5*2 near the SIT. Physically, the coupling to 
the diffusive degrees of freedom in the 2DEG quenches 
charge-fluctuations in the IDJJA, resulting in an en- 
hanced effective charging energy Eq. This phenomenon 
is more precisely described by the RG-flow equation for 
K (1) [cf. Eq. ([21])], where K (1) is renormalized to lower 
values by r] (1). Indeed, near the SIT, a small increase in 
the initial value r]o = r] {I = 0) (i.e., an increase in R\j) 
can effectively control the RG-flow of K (l) and therefore, 
that of A(/), inducing the SIT. We illustrate this point 
in Fig [2j where the schematic phase diagram obtained 
by integration of the RG-flow Eqs. (21)-(24), with initial 
parameter Aq = A (/ = 0) = 0.01. Note the stabilization 



6 



of the insulating groundstate due to Ohmic dissipation 
induced by the coupling to the 2DEG. 

In a first approximation, this effect is similar to the 
dissipation-driven SIT obser ved in 2DJJAs capacitively 
coupled to a diffusive 2DEG.^ ^ * ^^ * ^^ * However, important 
differences appear with respect to the 2D case. In that 
case, it was argued that dissipation produced a renormal- 
ization of the effective parameters Ej and Ec of the array 
due to the incomplete screening of the Coulomb interac- 
tion in a certain frequency-regime.^ Physically, the slow 
diffusive response of the 2DEG cannot follow the faster 
dynamics of the 2D plasma mode, and cannot screen it 
efficiently. However, in the ID geometry the ID plas- 
mon is effectively very well screened by the 2DEGPS1, and 
it could be naively concluded that no dissipation-driven 
SIT should be observed. However, this screening effect 
is compensated by the presence of strong backscattering 
occuring in ID [i.e., action 6*2 , Eq. (20)], and originated 
in the retarded interaction I'eff t). The net result is 
that the dissipation-driven SIT is restored in ID. 

Although one expects the nature of the non- 
superconducting groundstate to be of the Mott-insulating 
type, by analogy with the well-known results for the 
sine-Gordon model^ strictly speaking we cannot 

extrapolate the results in this Section to the strong- 
coupling situation, and a different method is needed in 
that regime. 



B. Self-consistent harmonic approximation 

To gain more insight into the phase in which the pa- 
rameter A (/) flows to strong-coupling, in this Section 
we make use of the variational self-consistent harmonic 
approximation^^. This method consist in finding the op- 
timal propagator g^^^ (q) of a Gaussian trial action of the 
IDJJA 



5'tr = ^Xl^tr^ (q) l^ql^ 



(25) 



where is the Fourier transform of ^(x,r). Here we 
have introduced the compact notation q = (/c, —uJm/u). 
The idea is to minimize the variational free-energy 



^var = ^tr + (5'eff " 5'tr)tr •> 

where the "trial" free-energy Ftr is 

Ftr = -^Y.^og[pLg,, (q)] 



(26) 



(27) 



The factor 1/2 in Eqs. ([25| and ([27) come from the 
constraint (j)* (q) = (j) (— q) since (j) {x~t) is a real field, a 
fact that reduces the number of independent degrees of 
freedom . 




Figure 3: Dimensionless SCHA parameters A and as a func- 
tion of calculated for the parameters 77 = 0.01, A = 0.05. 
We obtain non-vanishing values only in the region K < 2 
(note the abrupt increase in that region), consistent with the 
results of the RG-analysis. 



The minimization of Fvar in Eq. ( |26| ) with respect to 

): 

9tr^ (q) 



9tT (q) yields the self-consistent equation for ^tr (q) 
1 9 u ^> 



TTUK 



ttK 



St] cos {uJmT) 

' ar — 



aro 
- 1 



X e 



(28) 



In general, the solution of this equation has to be found 
numerically. However, for small A and rj the analytical 
solution 



Qtr^ (q) = 9ll (q) 

is obtained. Here gj^j^ (q) 



c 



A 

OTo ' 



(29) 



tK 



k'^ is the Lut- 

tinger liquid propagator, corresponding to the action Eq. 
(13) (for 7^ = 0). Physically, this propagator describes 
an insulator (given by a non-vanishing gap or "mass" 
term A) with Ohmic-dissipative dynamics (encoded in a 
non- vanishing Q. Note that dissipation dominantes for 
frequencies \uOm\ > A/Ctq. These parameters are found 
solving the following set of non-linear equations 



C Sttt] 



2K 



■ 



A = 4A 



(^K-K + 2Vk 



K 



(30) 
(31) 



obtained replacing the solution Eq. (29) back into Eq. 



28). Starting from the self-consistent solution of Eqs. 

30) and (31) for A in absence of dissipation (i.e. 7^ = 0), 
we can study the regime 77 <C A <C 1 perturbatively in 77, 
and we obtain the following estimate for the gap increase 
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due to dissipative effects 



SA ~ 27r' 



2VKAI 



(32) 



This result is consistent with the fact that dissipation 
in the density (i.e., field (j)) quenches charge-fluctuations 
and therefore favors an insulating groundstate. 

In Fig. [3] we show numerical results for A and ( as 
a function of K for the values A = 0.05 and r] = 0.01. 
Note the sharp increase of both A and ( for K < 2. This 
result is consistent with the RG-analysis, which predict 
the breakdown of the Luttinger liquid phase for K < 2 
in the weak-coupling regime. Within the SCHA, the 
physics of the strong-coupling fixed point is encoded in 
non- vanishing values of ( and A, providing a complemen- 
tary description to the RG-analysis. 



IV. TRANSPORT PROPERTIES 

In this section we concentrate on the dc-conductivity 
of the IDJJA, a quantity of central interest in 
experiments.^ ^ We first focus on the current-density 
j (x). Since the field \/0(x)/7t is the momentum of 
Cooper-pairs [cf. Eq. Q], the usual minimal cou- 
pling procedure {x) /tt ^ [V^ {x) — 2eA {x)] /tt (with 
e the electron charge and A the vector potential) in 
Hamiltonian Eq. ([2| allows to obtain the current as 
j (x) = —dHjjA/SA{x). In our problem, it explicitly 
read^ 

j {x) = uK [Ve {x) - 2eA (x)] . (33) 

The conductivity along the wire is obtained from the 
Kubo formula^^ 



cr(cj) 



xg(Q,c^) 

z (cj + i5) ' 



where we have restored the Planck constant and where 
we have used that y^-- (q) ^ in the limit k = 0.1^ Note 
that the real part of ctq (cj) consists of a Drude-peak at 
cj = 0, as expected for a superconductor. This result 
can be understood from the fact that the total charge 
current Je = ^ dx j (x) is a conserved quantity in absence 
of QPS and dissipation processes, i.e., it commutes with 
the hamiltonian Hjja- 

The effect of a finite rj in the Gaussian sector of the 
theory [cf. Eq. (13)] has been studied in Ref. 25, and 
produces a broadening of the plasmon peak, whose width 
r vanishes as T ~ \k\. Consequently, only taking into 
account this effect, a well-defined Drude-peak in a (uj) 
for cj = is recovered, and the system should behave as 
a perfect conductor. 

Let us now study the effects of the terms and 82- 
When A and rj are irrelevant perturbations (in the RG 
sense), their effects on the conductivity can be studied 
within the theoretical framework of the memory function 
formalism.HS In this approach, the central assumption is 
that the Kubo formula for the conductivity Eq. ( 34 ) can 
be recasted a^^ 



i{2ey 



uK 



■M{uj,T)' 



(36) 



where M (cj, T) (i.e., the memory function) is a meromor- 
phic function depending on the terms in the Hamiltonian 
responsible for degrading the current, and hence produc- 
ing a finite resistivity. Current-decay originated in QPS 
and in the coupling to the dissipative modes in the 2DEG 
induce finite resistivity in the IDJJA for all temperatures 
T < Tc. In particular for temperatures T <^Tc, and per- 
turbatively in A and 7^, we obtain 



0iT) = 



a{2ey 



[AiT 



2K-3 



(37) 



(34) 

where 



where xfj{k,^) = liiRi^^^^^is XjA^) is the re- 
tarded current-current correlation function and Xjj (q) = 
(i*(q)j(q)) = SHnZ/SA{q)SA*{q)\^^^ is the 
current-current correlation function obtained in the 
linear-response regime. It is convenient to express this 

-yd I 

^33 ^ ^33 



correlator as Xjj (q) = X% + X?^ (q)' where x'jj 



(2e) uK/n is the diamagnetic contribution and 



Ai = X^^TT^ 



K 



(38) 



A2 = 7/327r^ cos [(1 + i^) tt] 5 [1 + K, -1 - 2K] 



2Tia 



2K 



U J 

(39) 



{uKfe{e{c^)e{-ci)), (35) 

TT 



is the paramagnetic term.f^ In absence of current- 
decaying mechanisms [i.e., A = 77 = in Eq. ([6|], the 
conductivity writes 



S (cj) + iF 



where the function B {x^y) is defined SiS B (x^y) = 
r (x) r (y) /r (x + y\ and T {x) is the standard Euler's 
Gamma function.l^S The term - T^^"^ in Eq. (IstI) is the 



contribution due to QPS processes, consistent with for- 
mer theoretical predictions.^^ The second term ~ 
originates in backscattering effects induced by dissipa- 
tion, and is consistent with the behavior predicted by 
Cazalilla et alW^ This last effect can be interpreted as 
a frictional drag produced by the diffusive modes in the 
2DEG.^ Note that at lowest order in A and 77, the two 
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contributions add up independently, indic ating th at for 
temperatures T"" <T <^Tc, where T* = ^ Ai/ A2/2t:tq, 
the resistivity in the IDJJA is dominated by frictional 
drag, while for T < T* <^ Tc the effect of QPS takes 
over. 

The non-trivial effects due to the renormalization of 
the bare couplings can be taken into account integrating 
the RG-flow Eqs. ([2T| -(24), and injecting them in the 
above Eqs. ([38f and (|39|. We integrate the RG- 

fiow up to the scale given by the temperature a{l) = 
a (0) = u{l) /27rT, and we use formula Eq. (37) with 
the parameters of the model calculated at the scale a (1). 
This allows to obtain g (T (/)) vs T (/). 

In Fig. [I] we show the resistivity g{T) of the IDJJA, 
calculated for different values of the parameter K and 
using the estimations for the bare parameters Aq = 0.01 
and 7^0 = 0.01. The results are normalized to a "high- 
temperature" resistivity g (Tq), where Tq a (0) /ii = tq, 
represents a high-temperature cutoff in the theory (e.g.. 
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Figure 4: Dc-resistivity q (T) of the IDJJA, normalized to a 
"high-temperature" value ^(Tq), as a function of T/Tq, cal- 
culated for the parameters Ao = 0.01 and 770 = 0.01, and 
for different values oi K = tty^^j/^- A low-temperature 
upturn of g (T) signals the formation of the insulating phase. 

Note that for the values K = 2.5 and K = 2.3, the re- 
sistivity shows a monotonically decreasing behavior, in- 
dicating a superconducting groundstate and consistent 
with the RG-analysis of Sec. |III A We also note a small 
kink around T* ~ 0.4 Tq, signalling the aforementioned 
crossover from dissipation-dominated to QPS-dominated 
resistivity. For K = 2.1, the resistivity first decreases and 
then shows a low-temperature upturn, indicating that the 
array is near the quantum critical point Kc. Finally, for 
lower values of the insulating behavior in the IDJJA 
is clear. Since both the integration of the RG-fiow equa- 
tions and the calculation of the memory-function formu- 
las are perturbative in A and 77, the calculation of the re- 
sistivity must be stopped whenever A (/) or 77 (/) become 
of order unity. 



In Fig. [5j we show the resistivity as a function of 
T/To, calculated for fixed K = 2.3 and A = 0.01, and 
for different values of parameter rj. We see that for 
T] = (i.e., = in the 2DEG), the array shows su- 
perconducting behavior, and the resistivity due to QPS 
processes is well described by the predicted power-law 
g{T) - t2^-3^ with K = K {I ^ 00) 2.2 the renor- 
malized value predicted by Eq. ([23|. Upon increasing 
the parameter 77, the resistivity of the array increases, de- 
veloping the aforementioned kink, but most importantly, 
the low-temperature resistivity develops an upturn, indi- 
cating a dissipation-driven phase transition to the insu- 
lating phase. 




Figure 5: Resistivity of the IDJJA (in units of h/e^a ) as 
a function of T/To, calculated for parameters K = 2.3 and 
A = 0.01 and for different values of ry. Although in absence 
of dissipation, the array is in the superconducting phase, a 
dissipation-driven SIT occurs upon increasing ry, consistent 
with the results in Fig. [2] The curve g (T) ^ T^^"^ is shown 
for comparison. 

More insight into the insulating phase can be ob- 
tained using the Luther-Emery refermionization solution 
for K = l.^^S-^ In absence of dissipation (i.e., 7/ = 0) 
an exact solution is obtained in terms of non-interacting 
fermions, with a gap A = iraX in their spectrum of excita- 
tions. Using the Kubo formula, one obtains the following 
expression for the dc-conductivity at low temperatures 
T<C A 



Cr(cj) 



h 



2'nT _A/Tr/ \ 

— ' o{(jj). 



This contribution arises from the excited quasiparticles 
above the gap A, which have an exponential population 
at low enough temperatures. This infinite conductivity 
occurs because in absence of dissipation, excited quasi- 
particles are infinitely long-lived. Using the memory- 
function approach for the refermionized problem in the 
regime 7/ <C A, T <C A, we find the analytical result 
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a (w = 0) = 



2 



-A/T 



where C2 is a numerical coefficient C2 — As ex- 

pected, dissipation introduced a finite lifetime in the 
quasiparticles, and a finite resistivity is obtained at = 
0. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the properties of a linear JJA 
capacitively coupled to a diffusive 2DEG placed in close 
proximity. Using a bosonization approach, we have de- 
rived an effective model for the IDJJA, and have ob- 
tained its critical properties and phases at T = 0. Our 
main result is the possibility to observe a SIT tuned by 
the parameter r] ~ Rd/Rq [cf. Eq. (16)]. This setup 
could be used to investigate the superconductor-insulator 
transition in a IDJJA under better controlled experimen- 
tal conditions as compared to other setups used in the 



past. Our work could shed some light on the under- 
standing of other ID superconducting systems showing 
a similar behavior, such as ultra-thin superconducting 
wires built by molecular templating^ ^ or by e-beam 
lithography^ techniques. 

We have shown that besides the more or less trivial 
static screening effect, the presence of a 2DEG induces 
dissipative effects in the quantum dynamics of the IDJJA 
due to backscattering processes induced by the dynam- 
ically screened Coulomb interaction, and explicitly de- 
pend on the sheet-resistance R\j of the 2DEG. These dy- 
namical effects play an important role in the quantum 
phase diagram of the IDJJA. This situation is differ- 
ent from previous approaches in higher dimensions. 
Indeed, in ID the plasmon mode is almost statically 
screened^, and this would lead to the naive conclusion 
that dynamical effects are not important. However, a 
more careful analysis shows that backscattering origi- 
nated in the dynamically screened Coulomb potential 
has the effect of restoring the SIT. In our system, these 
dynamical effects have important consequences for the 
critical properties of the array, and should be possible 
to observe them in dc-transport measurements. Phys- 
ically, the coupling to diffusive modes in the metal in- 
duces charging effects which are local in space (i.e., of 
the order of the lattice parameter a of the IDJJA) but 
which are non-local in time (i.e., Ohmic dissipation ef- 
fects), and tend to quench charge fluctuations, rendering 



superconductivity weaker. 

By the means of a weak-coupling RG-analysis and a 
variational approach, we predict a SIT driven by the 
presence of dissipation in the 2DEG. This SIT is of the 
BKT-type and mediated by unbinding of QPS/anti QPS 
pairs, like in the dissipationless case.^ Near the critical 
line the effects of QPS are stronger than those originated 
in dissipation and results in a SIT. This scenario is cor- 
roborated by a subsequent variational analysis of action 
Eq. (12), which suggests the formation of a gap A in the 



spectrum of excitations of the IDJJA [cf. Eq. (29)]. 

Our results suggest that dissipation renormalizes the 
QPS-rate to higher values and the ratio ^/Ej/Eq to 
lower values [cf. Eqs. (21)-(24)], rendering supercon- 
ductivity in the IDJJA weaker. Eventually, an increase 
of R[j [and therefore of 7^, in view of Eq. ( 16 )] , could drive 
the system into the insulating phase, as can be seen in 
Figs. |2]and[5l This ph enomenon is different to the case 
studied by Cazalilla et a/., where QPS processes were ab- 
sent, and it was dissipation itself that drove the quantum 
phase transition for the critical value Kc = 1/2.^ 

We have also studied the consequences on the 
temperature-dependent dc-resistivity of the array g (T) . 
We have shown that a non- vanishing R^ induces a rich 
behavior of g{T). In particular in the superconducting 
phase, where the IDJJA is in the Luttinger liquid uni- 
versality class, and the effects of QPS and dissipation 
are perturbative, the resistivity of the array g{T) fol- 
lows a power-law behavior g (T) = AiT""^ + A2T''^ , with 
exponents ui = 2K — 3 and V2 = [cf. Eq. (37)] 
generated by QPS and dissipation, respectively. There- 
fore, the results of this paper could be relevant in the in- 
terpretation of experimental results of transport through 
superconducting circuits subject to dissipative effects. In 
the insulating phase, the low-temperature dc-resistivity 
is expected to show thermally- activated behavior.^ In 
particular for = 1, the resulting model can be stud- 
ied analytically with a refermionization approach, and 
results in a resistivity g{T) ~ rfFe^l^ j Quite impor- 
tantly, note in this expression that the resistivity depends 
also implicitly on ri via a renormalization of the gap 77. 



Acknowledgments 

We acknowledge useful discussions with I. Pop, W. 
Guichard and F. Hekking. This work was supported by 
the Swiss National Foundation under MaNEP and divi- 
sion II. 



^ G. Schon and A. D. Zaikin, Physics Reports 198, 237 (2001). 

(1990), ISSN 0370-1573. ^ D. B. Haviland, Y. Liu, and A. M. Goldman, Phys. Rev. 

^ R. Fazio and H. van der Zant, Physics Reports 355, 235 Lett. 62, 2180 (1989). 



10 



^ N. Mason and A. Kapitulnik, Phys. Rev. Lett. 82, 5341 
(1999). 

^ N. Mason and A. Kapitulnik, Phys. Rev. B 65, 220505 
(2002). 

^ A. Bezryadin, C. N. Lau, and M. Tinkham, Nature 404, 
971 (2000). 

^ K. Y. Arutyunov, D. S. Golubev, and A. D. Zaikin, Phys. 

Rep. 464, 1 (2008). 
^ A. Bezryadin, J. Phys.: Condens. Matter 20, 043202 

(2008). 

^ H. S. J. van der Zant, W. J. Ehon, L. J. GeerHgs, and J. E. 

Mooij, Phys. Rev. B 54, 10081 (1996). 
^° A. J. Rimberg, T. R. Ho, g. Kurdak, J. Clarke, K. L. 
Campman, and A. C. Gossard, Phys. Rev. Lett. 78, 2632 
(1997). 

Y. Takahide, R. Yagi, A. Kanda, Y. Ootuka, and S.-i. 
Kobayashi, Phys. Rev. Lett. 85, 1974 (2000). 
E. Chow, P. Delsing, and D. B. Haviland, Phys. Rev. Lett. 
81, 204 (1998). 

W. Kuo and C. D. Chen, Phys. Rev. Lett. 87, 186804 
(2001). 

H. Miyazaki, Y. Takahide, A. Kanda, and Y. Ootuka, 
Phys. Rev. Lett. 89, 197001 (2002). 

Y. Takahide, H. Miyazaki, and Y. Ootuka, Phys. Rev. B 
73, 224503 (2006). 

M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 
M. P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 
A. D. Zaikin, D. S. Golubev, A. van Otterlo, and G. T. 
Zimanyi, Phys. Rev. Lett. 78, 1552 (1997). 
J. E. Mooij and C. J. P. M. Harmans, New Journal of 
Physics 7, 219 (2005). 
^° L M. Pop, L Protopopov, F. Lecocq, Z. Peng, B. Pan- 
netier, O. Buisson, and W. Guichard, Nature Physics 6, 
589 (2010). 

V. E. Manucharyan, N. A. Masluk, A. Kamal, J. Koch, 
L. L Glazman, and M. H. Devoret, cond-mat/ 1012. 1928 
(2010). 

22 N. Giordano, Physica B 203, 460 (1994). 
22 A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 211 
(1981). 

2^ K.-H. Wagenblast, A. van Otterlo, G. Schon, and G. T. 

Zimanyi, Phys. Rev. Lett. 79, 2730 (1997). 
2^ A. M. Lobos and T. Giamarchi, Phys. Rev. B 82, 104517 

(2010). 

2^ M. A. Cazalilla, F. Sols, and F. Guinea, Phys. Rev. Lett. 
97, 076401 (2006). 



2^ P. Goswami and S. Chakravarty, Phys. Rev. B 73, 094516 
(2006). 

2^ G. Refael, E. Dernier, Y. Oreg, and D. S. Fisher, Phys. 
Rev. B 75, 014522 (2007). 

2^ T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, Oxford, 2004). 

2° V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971). 

2^ J. M. Kosterlitz and D. J. D J Thouless, J. Phys. C: Solid 
State Phys. 6, 1181 (1973). 

22 M. Tinkham, Introduction to Superconductivity (McGraw- 
Hill, 1996), 2nd ed. 

22 J. Negele and H. Orland, Quantum Many Particle systems^ 
Frontiers in Physics (Addison- Wesley, Reading, Mass., 
1987). 

2^ G. D. Mahan, Many-Particle Physics^ Physics of Solids and 
Liquids (Kluwer Academic/Plenum Publishers, New York, 
2000), 3rd ed. 

2^ E. Akkermans and G. Montambaux, Mesoscopic Physics of 
Electrons and Photons (Cambridge University Press, Cam- 
bridge, 2007). 

2^ M. Abramowitz and L Stegun, Handbook of mathematical 

functions (Dover, New York, 1972). 
2'^ F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

28 R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 

29 R. M. Bradley and S. Doniach, Phys. Rev. B 30, 1138 
(1984). 

L. L Glazman and A. L Larkin, Phys. Rev. Lett. 79, 3736 
(1997). 

A. Vishwanath, J. E. Moore, and T. Senthil, Phys. Rev. B 
69, 054507 (2004). 
^2 R. p. Feynman, Statistical Mechanics (Benjamin, Reading, 
MA, 1972). 

^2 W. Gotze and P. Wolfle, Phys. Rev. B 6, 1226 (1972). 
T. Giamarchi, Phys. Rev. B 46, 342 (1992). 
A. G. Rojo, J. Phys.: Condens. Matter 11, R31 (1999). 
A. Luther and V. J. Emery, Phys. Rev. Lett. 33, 589 
(1974). 

^'^ C. N. Lau, N. Markovic, M. Bockrath, A. Bezryadin, and 
M. Tinkham, Phys. Rev. Lett. 87, 217003 (2001). 

^8 Y. Chen, S. D. Snyder, and A. M. Goldman, Phys. Rev. 
Lett. 103, 127002 (2009). 

^9 Indeed, in the static limit cOm = 0, the effective screen- 
ing potential writes Vscrik^O) = 2e^ Kq (2kd) / Cr and in 
the limit kd ^ compensates the logarithmic divergence 
of vi-D{k,0) [cf. Eq. ([9|], and yields limfc^o '^^eff (^, 0) = 
2e2 \n{2d/a) /e^^ 



